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Abstract. - We show results of simulations of a weakly driven four dimensional Edwards- 
Anderson spin glass, which present clear signatures of dynamical ultrametricity at low tempera- 
tures. The presence of a hierarchical organization of time scales is evident from the appearance 
of a triangular relation between correlations at three characterisitc long times in the asymp- 
totic limit of small drive, and also in the scaling form of stationary correlations in this regime. 
Recent results in the three dimensional case have been inconclusive showing the importance of 
dimensionality in tackling this delicate problem. 



Introduction. - Ultrametricity between points or states in the phase space of complex 
systems is a remarkable property which signals a hierarchical organization of the states Jj],|^]. 
In an N dimensional space, a square distance between two states a and /3 with coordinates 
sf = ±1, i = 1 . . . N, can be defined as: 

1 N 

i=l 

= l-q af} (2) 

where q a/3 = (1/N) ^ sfsf is the overlap between the states. If a, (3 and 7 are three states, 
then ultrametricity means that: 

d ai <max{d a l 3 ,d 13 ' 1 } (3) 

or cquivalently: 

q a i > 77im{g Q/3 ,/ 7 } (4) 

Ultrametricity is a distinctive feature of the phase space of the Sherrington-Kirkpatrick spin 
glass The search for ultrametricity in more realistic short range models has been acom- 
plished by searching for ground states of small systems H or thermalized states at not too 
low temperatures Q. 
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Signatures of ultrametricity can be found also in the dynamical behaviour of the system. 
An aging system may present a kind of dynamical ultrametricity ■ In this case correlations 
between states at three long times t\ 3> t 2 S> £3 should obey, asymptotically, the triangular 
relation: 

C(tx,t 3 ) = min{C(t u t2),C(t2,t 3 )} (5) 

It is difficult to test this relation in simulations due to strong preasymptotic time effects. 
An indirect test by constraining two real replicas to have one of the equilibrium values of 
the overlap was implemented in J9| for the three dimensional Edwards- Anderson spin glass. 
Another approach has been proposed in JlO| | by considering how dinamical ultrametricity 
should look like in a system subject to a weak external drive. When a system is subject to 
an external driving force, aging is eventually interrupted and stationary dynamics is attained 
after a characteristic time r £ , which depends on the strength of the drive, e. Interestingly, in 
the case when the three times are strongly separated, the relation (JsJ> can be translated to the 
stationary dynamics: 

C(ti - f 3 ) = C(ti - t 2 + 1 2 - t 3 ) = min{C(h - t 2 ),C(t 2 - t 3 )} (6) 

In the preasymptotic region one expects that, approximately: 

C{ti+h) = f[C{t x ),C{t 2 )\ (7) 

This approach is interesting because it is conceptually simple and, in principle, can be tested 
straitforwardly in computer simulations. In fact, in ref. jlO]], the method was applied to a 
driven Edwards- Anderson spin glass in three dimensions, with inconclusive results. In this 
letter we perform a similar study on the four dimensional EA spin glass with asymmetric 
couplings. Our results are very different from those obtained in three dimensions, and the 
presence of dynamical ultrametricity is clearly established. Our results suggest that, in order 
to detect ultrametricity in finite dimensional models, the drive must in fact be very small and 
consequently very long time scales need to be reached. This is at variance with the observed 
behaviour in mean field models. 

Four dimensional Edwards- Anderson spin glass. - The model is defined by the Hamilto- 
nian: 



H 



N 
<i,j> 



SiJijSj 



(8) 



where {si = ±l,i = 1 . . . iV} are N Ising spins and < i,j > denotes a sum over nearest 
neighbors. A drive can be applied to the system by adding a non Hamiltonian contribution to 
the energy. In a spin glass this can be accomplished, for example, by adding a non symmetric 
part to the copulings matrix. In this case the couplings Jij can be chosen as a weighted sum 
of a symmetric and a completely asymmetric part pT| : 



J' i j > 

Vl - 2e- 



2e 2 



(1 



The symmetric part of the interaction is given by J^p — J^ J — ±1 with probability 0.5. The 



el. 



AS) 



(NS) 



(9) 



non-symmetric part is chosen independently of J^ D> . Finally, e measures the strength 

of the non-symmetric part (the drive). 



(NS) 
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Fig. 1 - Stationary autocorrelations for three values of the drive, from top to bottom e = 0.1, 0.15 
and 0.2. 



We let the system evolve until it reaches the stationary regime and then measured auto- 
correlations: 

1 N 

C(t 1 -t 2 ) = -J2si(ti)s i (t 2 ) (10) 
i=i 

The measures were done on a system of linear size L — 10 at temperature T = 1, approximately 
0.5T C , for three different values of the driving force e = 0.1, 0.15 and 0.2. 

In the limit e — > the correlations develop a plateau at a value C rs qEA- The development 
of the plateau can be observed only in the limit of very small drive, when the time scales for 
structural relaxation become very large, as can be seen in fig. [j]. This growth of relaxation 
time when e — > is equivalent to what happens in an aging system when the waiting time 
t w — > oo. 

Dynamical ultrametricity implies the existence of a hierarchy of time scales. In particular 
this implies that the correlation will stay at nearly constant values for successively longer 
times as the system relaxes. This should manifest ideally in the development of a series of 
plateaus in the correlation as discussed in fllPf . In such a scenario it is clear that a scaling 
of the autocorrelations assuming a single relevant time scale, r(e), should fail. This simple 
scaling means that 

c,W=</(4) (ID 

In fig. H we show that in fact this does not work in the four dimensional spin glass. We 
tried to scale the times with the time at which the correlations relaxed at C=0.5, this value 
being well below the plateau at qEA ~ 0.7. The three curves corresponding to the three 
different drivings meet only at C=0.5, by construction, and cross each other at that point. 
This behaviour is similar to what happens in systems with full replica symmetry breaking, 
like the Sherrington-Kirkpatrick model or the Hofield neural network model |^,|l2|, which are 
ultrametric. A stronger evidence of dynamical ultrametricity is the presence of a logarithmic 
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Fig. 2 - Simple scaling of the stationary autocorrelations for the three values of the driving force 
e — 01, 0.15 and 0.2. Times are rescaled by the values at which the respective correlations decay to 
0.5. 



scaling of the form 



C 6 (t) cx / 



ln(t) 



(12) 



>(r(e)), 

A scaling of this type is shown in fig. |[ It is much better than the previous simple scaling 
and is excellent for the two smaller drivings. This is an indication that the relevant results 
can only be reached approaching effectively the asymptotic regime of e — > 0. 




log(t)/log(t(0.5)) 



Fig. 3 - Logarithmic scaling of the stationary autocorrelations for the three values of the driving force 
e = 01, 0.15 and 0.2. 
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It is easy to show that a scaling of the form (|i2|), in the presence of a small driving force, 
implies dynamic ultrametricity [[To). In fact, eq.([12|) implies 

t(C e ) ex (T(e)f C ^ (13) 

where j(C e ) is the inverse of /, a positive decreasing function of the correlation. Considering 
two values of the correlations such that C\ < C2 < q_EA one finds that: 

^1 « (r(e))^-^) ==S 00 (14) 

When, as in this case, the two times become assymptotically infinitely separated, this 
together with eq.(^) implies that C(t(C\)) — min{C(t(Ci)),C(t(C2))}, which is a realization 
of the ultrametric relation eq. (^|) as applied to the stationary dynamics at finite e. 

We have performed an analysis similar to the previous one to the four dimensional EA spin 
glass with gaussian couplings and without drive, i.e. at e = 0. The resulting aging dynamics 
in that case was studied in . The two time autocorrelation was found to scale as: 

ffh=*a\ 

C(h,t 2 ) = {qea + a(h - t 2 )~ x ) K f 2 1 (15) 



with the scaling fuction: 



1 constant for z — > 

/(Z) H Z ~A(T) for z _ QO _ (16) 



This scaling gives an autocorrelation in the aging regime that is not ultrametric, at variance 
with our results in the driven system. The difficulty in observing ultrametricity during aging 
dynamics was already noted in ||. Nevertheless, Eq.(|l5|) is not the only scaling form which 
correctly fits the data for the autocorrelation functions in the aging regime. Also ultrametric 
scaling forms can be used with equivalent results |l4| ]. Clearly the definitive scaling during 
aging is still not known. 

As a last test, we have calculated the function f(C\,C2) of eq.(^). As we are proving 
the preasymptotic regime, this function can give a flavour of the tendency of the triangular 
relation to become ultrametric or not. In this driven dynamics scenario, ultrametricity means 
that, as e — > 0, f(C\, C2) — > min(C\, C2). In other words, the curves of constant / projected 
in the (Ci, C2) plane should tend to have right angles 10 1. In fig. ^ we show our results for 
the 4D EA spin glass. 

Each group of three curves correspond to a fixed value of / and the three values of the 
drive parameter studied. In each of these groups upper lines correspond to e = 0.2, middle 
lines to e = 0.15 and lower lines to e = 0.1. The last group shows only two lines corresponding 
to e = 0.2 and e = 0.15 because of the lack of data for the smaller e = 0.1 in this region 
of very small correlations. This figure should be compared to the corresponding one for the 
three dimensional Edwards- Anderson model of reference |Io(| . In the four dimensional case 
the separation of the lines becomes stronger on longer times. Also, for the lower values of / 
a well defined tendency towards right angles is seen upon decreasing e, i.e. correlations are 
flowing in such a way as to satisfy dynamical ultrametricity as the drive e — > at long times. 

In the framework of the rheological approach to glassy dynamics, we have presented strong 
evidence of an ultrametric structure underlying long time correlations in the four dimensional 
+-J Edwards- Anderson spin glass. In comparing our results with the corresponding ones on 
the three dimensional model of |l(| some comments are in order: first, we tried to get closer 
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Fig. 4 - Projections on the plane (CijCa) of some constant values of the function / defining the 
triangular relation eq.( (jj), for three different values of the drive intensity e. Each group of three 
curves correspond, from top to bottom, to e = 0.2, 0.15 and 0.1 (see text) 



to the asymptotic regime e — > in which analytical results are valid. Note that, in fact, the 
logarithmic scaling of fig. |^ is very good only for e < 0.2. Of course, the price to pay is a 
strong growth of the relaxation time on decreasing e, as can be seen in fig. 0. Also, m mean 
field models relaxation is faster and relevant results can be seen already in the preasymptotic 
regime, i.e. for not too small values of e. In three dimensions relaxation is still slower than 
in four dimensions because of the larger value of the dinamical exponent z. To have an idea 
of the time scales involved note that z(T c ) = 4 in the mean field model, z(T c ) = 5.26 in four 
dimensions and z(T c ) — 6.67 in three dimensions. The difference is roughly 25% between 3 
and 4 dimensions and 30% between 4 dimensions and mean field. The relaxation times in 
three dimensions grow with respect to four dimensions as: 

T 3d (T/T c ) » T id (T/T c y 3d/zid (17) 

While at present it is a very hard problem, the three dimensional case should be accesible 
in the near future. 

Once more, the four dimensional spin glass turns out to be a very valuable tool for ad- 
dressing the question as to what extent some characteristics of the rich phenomenology of 
mean field spin glasses are present in finite dimensions below the upper critical one. Finally, 
probably results similar to ours can be observed in three dimensions on much longer time 
scales than those studied up to now. 
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